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1. Introduction 



Abstract. The dependence of the strangeness satm^ation factor on the system 
^ , size, centrahty and energy is studied in relativistic heavy-ion colhsions. 
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Statistical-thermal models (see for a review) aim at reproducing hadron multi- 
Oh! plicities in relativistic heavy-ion collisions in certain phase space regions by a small 

number of parameters. As these parameters evolve smoothly with external condi- 
tions, thermal models have predictive power. Moreover, the thermal parameters, 
taken literally, allow for a quantification of the conditions reached in the transient 
stage of compressed and heated strongly interacting matter. While there is exhaus- 
tive literature on the successful use of the thermal models, the foundation of the 
I applicability is far from settled. One should also keep in mind that thermal models 

address only a subset of the wealth of observables. 

The present contribution summarizes our recent work El 0] investigating the 
dependence of strangeness saturation on system size, centrality and energy, both in 
full phase space (47r) and at mid-rapidity (y w 0). 

2. Thermal Model 

Our model is based on the grand canonical ensemble expression for the abundance 
of hadron species i (Fermions [+1] or Bosons [-1]) with degeneracy gi 
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where T is the temperature, /* a set of chemical potentials corresponding to the 
conserved charges Q; Ei — -^/p^ + mf with vacuum mass m.^. 

For small particle numbers one has to turn to a canonical ensemble, e.g. by 
the projection method It is important to include resonances (we include me- 
son (baryon) states composed of the three lightest quarks and anti-quarks up to 
2.3 (2.6) GeV); their widths are taken into account via the Breit-Wigner param- 
eterization BW{mi,Vi), which collapses to a (5-function for the respective ground 
states. Due to feeding and resonance decays, the final hadron multiplicities become 
Ni = Nf'"" -f \Nf'"'Br{j ^i)- Nf"''Br{i where Br denotes the 

corresponding branching ratio. 

Our focus here is the phenomenological strangeness saturation factor 7^ in 
Eq. (1). It is thought ['S' to parameterize possible deviations from chemical equilib- 
rium in the strange sector. Whether one common factor 7^ is sufficient, with powers 
\Si\ determined by the total strangeness content of hadron i, needs still to be tested 
by analyzing data. Also in the non-strange sector such an off-equilibrium parameter 
should be introduced [^_. Due to the small data samples at our disposal, such an 
additional parameter cannot be unambiguously fixed; therefore, we are forced here 
to assume full equilibration in the non-strange sector. 

The parameters T and appear exponentially in Eq. (1), while the sensitivity 
on variations of 7^ is weaker. For this reason, we study the mentioned parameter 
dependencies based on analyses of data sets comprised of the same hadron species. 
Otherwise variations of 7^ can be absorbed in small changes of T and fi. Unfortu- 
nately, this shrinks the wealth of available data to quite poor subsets. 

Within the Cooper-Frye formalism, dynamical effects factor out of fully- integrated 
hadron yields [IS], provided that freeze-out occurs on a hyper-surface with constant 
T and fl. In this respect, Eq. (1) covers also dynamical situations. Instead of consid- 
ering hadron ratios, we consider multiplicities A^f"'™ with a fiducial normalization 
volume V in Eq. (1). 

3. Data Analysis 

3.1. CERN-SPS energies 

3.1.1. Centrality and System Size Dependence at -Ebeam = 158 AGeV 

A-K multiplicities of tt^, , cf) and A^part are at our disposal for (i) central collisions 
of C + C and Si + Si, and (ii) centrality binned collisions Pb -I- Pb. More baryon 
information is desirable, but only p's are available additionaly for Pb + Pb. The 
left panel of Fig. 1 summarizes our findings [ 3: . We observe, 

• Pb + Pb: 7s increases with centrality, but stays below unity, irrespective of 
whether p's are included or not; inclusion of p's reduces 7^ somewhat; minimizing 

or the quadratic deviation of data and the model multiplicities, g^, yields slightly 
different results. 

• C 4- C, Si + Si: the extracted values of 7s off-set the down-extrapolation for Pb 
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+ Pb. This again emphasizes the findings of NA49 [ 9;: system-size and centrahty 
dependencies differ. 

• Most striking is that 7s goes as the number of participants which underwent 
multiple collisions, as extracted from a Glauber model calculation (cf. Fig. 3 in [ 

El)- 

There is some controversy on the use of 47r or mid-rapidity data. Our analysis [ 
0] of TT^ , and centrality-binned mid-rapidity data from Pb + Pb at 158 AGeV 
reveals that 7s does indeed differ for An and y « data (see Fig. 1, middle panel). 
This is due to the fact that the rapidity distributions of various hadrons differ. 
Combining the results, even the very restricted data samples point to 7s = 1 ± 0.25 
for most central collisions. This suggest that 7s = 1 with slightly renormalized 
values of T and /2 is a good choice for the given beam energy and high centrality. 

3.1.2. Energy Dependence 

Unfortunately, the previous analysis cannot be extended to extract energy depen- 
dencies. Rather, the An multiplicities of ir^, K^, A and A at beam energies of 
40, 80 and 158 AGeV are at our disposal as comparable data samples. The re- 
sults of our analysis [01 are shown in Fig. 2: 7s displays a dip at 80 AGeV, while 
7s (40 AGeV) - 7s (158 AGeV). (The significance of the drop at 80 AGeV stiU needs 
to be clarified.) 7s(158AGeV) is compatible with the strangeness saturation factor 
reported in the previous subsection emerging from a data set suitable for a central- 
ity analysis. The open circle follows from an analysis of mid-rapidity data of tt^, 
K^, A and A. The change relative to results reported in the previous subsection 
highlights the sensitivity on different data sets, if the number of analyzed hadron 
species is small. 

3.2. RHIC Energies 

Our analysis is restricted to the published data of tt^, and at y w in 
centrality-binned coUisions of Au -I- Au at y^JJvjv — 130 GeV The upcoming 

data for further hadron species, including excited states, and data at y^saT/v = 200 
GeV will complete our picture on the systematics of hadron multiplicities. 

3.2.1. Centrality Dependence 

As for CERN-SPS energies, 7s increases with centrality, see right panel in Fig. 1. 
The overall trend follows again the number of multiple participant collisions ob- 
tained in a Glauber model calculation The most central collisions point to 
strangeness saturation, in particular when feed-down from weak decays is included. 
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3.2.2. Energy Dependence 

Fig. 2 suggests that increasing energy causes increased strangeness saturation. How- 
ever, this conclusion should be taken with caution, as Fig. 2 displays in one plot 
analyses of An and j/ « data and data sets comprising different hadron species. As- 
suming that at CERN-SPS the An data are sensible, while at RHIC only y « data 
should be analyzed within thermal models, then from Fig. 2 one does indeed derive 
the mentioned expectation: with increasing energy, 7^ evolves towards saturation. 
One observes in Fig. 2 that the thermal model delivers for pp and pp collisions 
values of 7s compatible with peripheral collisions. 
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Fig. 1. Left panel: System-size and centrality dependencies of 7s, as extracted 
from centrality-binned Pb-|-Pb [El E| and central C-l-C and Si-|-Si data [E] under 
various fit conditions, assuming 50% feeding from weak decays. Also shown is the 
fraction of participants which underwent multiple collisions, /2. Middle panel: 
Comparison of 7s extracted from mid-rapidity NA49 data [ with the results of 
our earlier analysis of NA49 47r-yields Right panel: 7s observed in Au-|-Au 

collisions as extracted from PHENIX data [Eli /2 from our Monte Carlo based 
Glauber model calculation [i3i. 



4. Agreement with Data 

Up to now we have reported the system-size, centrality and energy dependence of 7s. 
The values of T coincide approximately with the chiral (deconfinement) transition 
temperature of T ~ 165 MeV, smoothly increasing from 150 MeV at beam energy of 
40 AGeV. fiB (the baryon potential dominating fl) decreases from 350 to 220 MeV 
at CERN-SPS and drops to about 30 MeV at RHIC. The actual values depend 
on the data samples analyzed. The agreement of data with multiplicities from the 
thermal model is extremely good with the following exceptions: 
• (j)/K~^ is underpredicted for peripheral collisions at beam energy of 158 AGeV, 
even if the (p yields are included in the fit; 
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• A + S° yields follow the data trend but are in absolute normalization below 
data when using the thermal parameters extracted from tt^, and p^ yields at 
y/sNN = 130 GeV; the S~ yields are at the lower end of the error bars. 

To understand which trends in the data drive the change of 7^ we exhibit in 
Fig. 3 scaled multiphcities as a function of -/Vpart. A scaling of p by A^pait results in a 
flat distribution at beam energy of 158 AGeV. The scaled tt* is slightly decreasing, 
while increases. It is the relative decrease of the "gap" between pion and 
kaon yields which drives 7^ towards unity. Similarly, at ^/snn = 130 GeV, the p^ 
distributions are flat, as also K^, but the tt^ decrease with 

increasing A^p art • Again, 

the decreasing "gap" between tt^ and drives the increase of 7^. Note that the 
effects are small and hardly visible in a log plot with unsealed yields. 

As an aside, we mention that we adjusted the parameters of the color glass 
model [ and obtain a good description of tt^ and p^ yields, while the model 
deviates in details from the data, see right panel of Fig. 3. 

5. Summary 

In summary we report here the dependence of the strangeness saturation factor, 
7s, on system size, centrality and energy. Most striking is the increase of 7s with 
centrality and the number of multiple participant collisions. Details of the energy 
dependence of 7s are puzzling and hampered by restricted data sets. Whenever 
possible, one should compare the findings with results obtained from the richer 
data samples for central collisions. Clearly, the present investigations should be 
extended to include multi-strange baryons and resonances in the analysis of more 
complete data sets. This will further test the applicability and predictive power of 
thermal models of statistical hadronization in describing hadron multiplicities. 
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Fig. 3. Scaled hadron yields as a function of A^part- Left (right) panel: Pb + 
Pb (Au + Au), symbols: data El Ej, dotted and solid lines: results of the 
thermal model, dot-dashed curves: our results for the color glass model [ITS). 
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